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Abstract 



We study crossover phenomena in a model of self-avoiding walks with 
medium-range jumps, that corresponds to the limit N — > of an A-vector 
spin system with medium-range interactions. 

In particular, we consider the critical crossover limit that interpolates be- 
tween the Gaussian and the Wilson-Fisher fixed point. The corresponding 
crossover functions are computed using field-theoretical methods and an ap- 
propriate mean- field expansion. 

The critical crossover limit is accurately studied by numerical Monte Carlo 
simulations, which are much more efficient for walk models than for spin 
systems. Monte Carlo data are compared with the field-theoretical predictions 
concerning the critical crossover functions, finding a good agreement. We 
also verify the predictions for the scaling behavior of the leading nonuniversal 
corrections. 

We determine phenomenological parametrizations that are exact in the 
critical crossover limit, have the correct scaling behavior for the leading cor- 
rection, and describe the nonuniversal crossover behavior of our data for any 
finite range. 



PACS Numbers: 05.70.Jk, 64.60.Fr, 75.10.Hk, 61.25.Hq 
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I. INTRODUCTION 



The universality of critical phase transitions is related to the presence of a diverging 
correlation length £. When £ is much larger than any microscopic scale characterizing the 
system, one observes a scaling behavior that is universal, i.e. independent of the microscopic 
details. However, in experimental situations the correlation length may not be so large and, 
on the contrary, it may be comparable to some other scale intrinsic to the system. In this 
case, one does not observe the expected critical behavior, but rather a crossover. Here, 
we will be interested in the crossover between the standard Wilson-Fisher behavior (near 
the critical point) and the mean-field behavior (far from the critical point) that is observed 
by varying the temperature in systems belonging to the iV-vector universality class, i.e. 
magnets, fluids, multicomponent fluid mixtures, ... (the critical behavior of these systems 
is reviewed, e.g., in Refs. Such a crossover is characterized by the Ginzburg number 

G || that measures the relevance of the magnetization (or density) space fluctuations that 
determine the departure from the Landau mean-field behavior. If t = {(3 C — /3)/P c is the 
reduced temperature, for \t\ > G the system shows an approximate mean-field behavior, 
while for \t\ < G one observes the standard Wilson-Fisher criticality. The crossover behavior 
is nonuniversal since it depends on the specific details of the system under investigation, 
and is usually described in terms of phenomenological models (see, e.g., Refs. Hj-jOJ and 
references therein for a discussion of phenomenological models for fluids, binary mixtures, 
and polymers). However, in a specific limit — we call it critical crossover limit — one can 
define universal quantities that do not depend on the microscopic details. 

In this paper, we will consider spin models with medium-range interactions. For instance, 
we may consider the Hamiltonian 

H = ~Y1 2 Vx'Vy, (1) 

x y.\y-x\<R 

where a x is an iV-component vector satisfying a x ■ a x = 1. The crossover behavior of 

these systems has been extensively studied numerically |JT~7| |22] , p~3[] . By means of scaling 



arguments, it was shown flT7|-|i9[| that the Ginzburg number G is proportional to _R~ M /( 4_d ) 
in d dimensions. Thus, for \t\ 3> R~ 2d /^ A ~ d ^ such systems show an approximate mean-field 
behavior, while for \t\ -C i?~ 2d /( 4 ~ d ) one observes the standard Wilson-Fisher criticality. 

Such a crossover can be described by using effective exponents. For instance, one can 
define an effective susceptibility exponent 7 e ff(i) (often called Kouvel-Fisher exponent p3| ) 
by 

*«(«.*> = — W^. CO 

XR{t) at 

where Xnif) is the susceptibility. By varying the temperature, the exponent 7 e fr(£, R) varies 
between 1, the mean-field value, and 7, the Wilson-Fisher value (actually, the full crossover 
behavior can be observed only for R large enough, in lattice models for R > 3). At R fixed the 
crossover behavior is not universal, and therefore the function 7 e ff(i, R) at fixed R cannot be 
predicted without an explicit reference to the microscopic details of the system. In this case, 
if one wishes to obtain interpolations that provide reasonably precise approximations, one 
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must resort to phenomenological models, such as those presented in Refs. p]-|l5|. However, 
there is a particular limit — the critical crossover limit — in which the effective exponents 
become universal. If we consider the limit t — > 0, R — > oo with t = t/G ~ £_R 2d /( 4_d ) fixed, 
then 7efr(t, -R) converges to a critical crossover exponent 7 e ff(t) that is universal apart from a 
trivial rescaling of t. In practice, at least in lattice models such as (|1|), the critical crossover 
functions provide a good description of the crossover as soon as the interactions extend over 
a few (two or three) lattice spacings. 

The universality of the critical crossover functions can be shown explicitly in the large- N 
limit ]2"3],|2"5| and for any value of N by performing an expansion around mean field 



Moreover, as shown in Ref. |[25| , they are given exactly by the field-theoretical crossover 
curves computed within the 4 framework in Refs. [pq-pO| . 

In this paper we wish to check with a high-precision simulation the field-theoretical 
predictions of Ref. p5 . Large-scale Monte Carlo results have already been obtained for 



the Ising model, both in two and in three dimensions (T^-f2T|. However, because of the 
difficulty of keeping under control the finite-size effects, only small ranges R were simulated 
near the Wilson-Fisher point. Although the general trend of the data was consistent with 
the analytic field-theoretical predictions, still in the Wilson-Fisher region a high-precision 
numerical check could not be done. The considered values of R were too small and there 
were significant discrepancies between numerical data and theoretical predictions. 

Here, we address the problem for a spin model in the limit iV — > 0, which can be described 
in terms of self-avoiding walks (SAWs) |3T]-[36|j. The advantage of such a system is that we 
can now work directly in the infinite-volume limit without finite-size effects and thus we 
can investigate systems with much larger values of the correlation length (in this paper 
we reach £ m 500 for systems in which the interaction extends up to 12 lattice spacings). 
In the limit iV — > 0, the model (|IJ) is mapped into a model of SAWs with medium-range 
jumps, i.e. of SAWs such that the length of each link is less than or equal to R. As usual 
in walk simulations, we work in a monodisperse ensemble, i.e. with walks of fixed length 
n. The length n replaces here the reduced temperature t. Medium-range SAWs show a 
crossover behavior depending on nG ~ nR~ 2d l <yA ~ d \ For n <C _R M /( 4_C the SAW behaves 
as an ordinary random walk (mean- field behavior), while in the opposite regime the self- 
repulsion becomes important and one observes the standard critical behavior. As we already 
stressed, for fixed values of R, such a behavior is not universal and can only be described 
phenomenologically. There is however a universal limit, the critical crossover limit: If we 
take the limit n — > oo, R — > oo keeping n = nR~ 2d l <yA ~ d ^ fixed, the crossover functions 
become universal, and can again be computed by using field-theory methods. 

The paper is organized as follows. 

In Sec. we introduce the model and the basic observables we consider. 
In Sec. [nj we discuss the critical crossover limit. In Sec. [HI A| we define the limit in spin 
models |[T8|-pO"]| and review the results of Ref. ||25|| . In Sec. [Ill B| we define the crossover limit 



for walk models and derive some general results for the universal crossover functions. In 
particular, we show that they are exactly related to the crossover functions computed in field 
theory p7|-^]. Moreover, by using the results of Ref. p5| , we show that by an appropriate 



definition of the range R the leading corrections to the universal crossover functions scale 
as R~ d , d being the dimension, as R — > oo. 

In Sec. [IV] we derive the expressions of the crossover functions from field theory general- 
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izing the results of Ref. |]38] . Details are reported in the Appendix, where we also compute 
the first coefficients of the asymptotic expansion of the crossover functions near the Wilson- 
Fisher point by using the fixed-dimension expansion in the zero-momentum scheme EDI and 



in the dimensional regularization scheme without e-expansion |^|29],|30| . 
In Sec. |V| we briefly describe the numerical algorithms we use. 

In Sec. |V]] we perform a detailed comparison of the numerical results in three dimensions 
with the field-theory predictions. We find a very good agreement, the deviations being small 
already when the interaction extends over three lattice spacings. Particular care has been 



devoted to the behavior of the leading corrections, 
predicted in Ref. p5 



We show that they scale as R as 



Finally, in Sec. [VII] we report our conclusions and discuss some further applications of 
these results. In particular, we give phenomenological expressions that are able to describe 
the crossover curves even outside the critical limit, for all ranges R we have considered. 

Preliminary results appeared in Ref. 



II. THE MODEL 

In this paper we consider SAWs with medium-range jumps. To be specific, let us consider 
a hypercubic lattice in d dimensions. Given an integer number p, let us define a lattice 
domain D p (x). If a; is a lattice point, D p (x) is the set of lattice points defined by 

D p (x) = jy : ~Vi\ < Pj ■ ( 3 ) 

We indicate with V p the number of points belonging to D p (x) and with R the mean square 
size of D p (x). Explicitly, we define 

V P = £ 1, (4) 



p yeD p (o) 



In three dimensions 



\/ p =I(2p + l)(2p 2 + 2p + 3), (6) 

2 = p(p + l) p 2 + p + 3 

10 2p 2 + 2p + 3' 1 ' 

For p — > oo, V p ~ |p 3 and R 2 ~ ^jp 2 . In the following we will often characterize the size of 
the jumps by R instead of p and thus we will write D R , Vr, . . ., instead of D p , V p , . . . 

Let us now define our model. We define an n-step _R-SAW as a sequence of lattice points 
{ou , . . . , u n } with ujq = (0, . . . , 0) and ujj + i G D R (ujj), such that Ui ^ Uj for all i ^ j. All 
walks are weighted equally. For p = 1 the model corresponds to a standard SAW with 
nearest-neighbor jumps. 
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We will consider the following observables: If c n ^(x) is the number of n-step i2-SAWs 
going from to x, we indicate with c nj # the total number of n-step walks and with E\ R the 
mean square end-to-end distance. They are defined as follows: 



c n,R 



^2cnA X )> ( 8 ) 



= — E^.^)- (9) 

Cn,R x 

This model of walks is related to a lattice A-vector model with medium-range interactions 
in the limit N — > 0. Indeed, consider the Hamiltonian 

x y<=D R (x) 

where a x is an A-dimensional vector satisfying a x ■ a x = 1, and define as usual 

Z^E^e-^, (11) 
M 

G R (x; (3) = (a ■ a x ) R = V a ■ a x e~^\ (12) 

ZdP) ^ 

The susceptibility and the (second- moment) correlation length are then defined as 

XK (/3)=^G R (x;/3), (13) 

X 

^25^5> 2cM ' (14) 

A standard procedure |3T|-|36|1 allows to prove that 

oo 

hm Xi? (/3) = ^/3 n Cri ,H, (15) 



N- 

n=0 



oo 



ra=0 

This equivalence will allow us to use the results available for the Hamiltonian flTOp that are 



discussed in detail in Ref. 125 



III. CRITICAL CROSSOVER LIMIT 

In this Section we derive some general results for the critical crossover limit of medium- 



range SAWs. They will be obtained by extending to walk models the results of Ref. [25|. 
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A. The variable-length /3-ensemble 



Let us consider the Hamiltonian fllOD , which, for R fixed, defines a generalized N- vector 
model with short-range interactions. For each value of R there is a critical point[] f3 c ,R, for 
P —> Pc,r the susceptibility and the correlation length have the standard behavior 



XR (P)^A x (R)t^(l + B x (R)t A + -- 
eM*As(R)t- 2 »(l + Bs(R)t A + -. 



(17) 
(18) 



where t = (j3 c ,R — P)/0 c ,r an d we have neglected additional subleading corrections. The 
exponents 7, u, and A do not depend on R. In two dimensions, for N = 0, 7 and v are 
known exactly 



3 

4' 



7 



43 
32 : 



(19) 



In three dimensions, for 



while A is still the object of an intense debate [[0,|7,|I),|£, 
N = 0, the best estimates of the exponents have been obtained in Monte Carlo simulations: 
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A 



0.5877 ±0.0006 
0.58758 ±0.00007 

1.1575 ±0.0006 



Ref. [50], 
Ref. P, 



Ref. [52 



(20) 



0.515; 



-0.017 
-0.007 



Ref. BT|. 



Less precise Monte Carlo results can be found in Refs. [^5| . [53| . |54l , ^6l and references therein. 
Similar, although less precise, results are obtained by using field-theory methods and from 
the analysis of enumeration series (for a list of results, see Refs. [p8| , |57| -|62]1 and references 
therein). 

On the other hand, the amplitudes are nonuniversal and depend on R. For R — > 00, 
they behave as 



A X (R) 
B X (R) 



gooj£2dA/(4-d) 



As(R) ~ A fR A{2 - dv V {A - d \ 
B^R) w BfR 2dA /^~ d \ 



(21) 



Corrections to these asymptotic behaviors vanish as |25| R~ d . 

The critical point (3 c ^r also depends on R. The expansion of /3 Cj r for R — ► 00 was derived 
in Ref. |25j in two and three dimensions. Explicitly, for = and d = 3, we have 



0, 



c,R 



1 

Vr 



l + I 



(22) 



where T\ and T2 are constants and Ir is a function of R. The nonperturbative constants t\ 
and T2 depend on the domain. The expression of T\ for a generic domain is reported in Ref. 



For the general iV-vector model in two dimensions, a critical point exists only for N < 2. Theories 
with N > 3 are asymptotically free and become critical only in the limit (3 — > 00. 
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TABLE I. Estimates of R 3 Ir for several values of p for the domain (|3|). From Ref. [25] 



p R 3 Ir p R 3 Ir 

3 0.043960387 10 0.043486698 

4 0.043921767 12 0.043451767 

5 0.043713672 14 0.043429899 

6 0.043664053 16 0.043415345 

7 0.043574469 18 0.043405187 

8 0.043547206 20 0.043397824 



2q| . For the domain ([!]), t\ ~ —0.00060(11). The constant t<i will be computed numerically 



in Sec. |VI A[ The function I R is defined by 



d*k 1 - U R (k) 

Ir = 1 wr^*r> (23) 



where 



u R (k) = i-l- elk ' x - ( 24 ) 

R xeD R (o) 

r R ~ crR~ 3 + 0(R~ 5 ). For the domain considered in this paper a ~ 0.04336529. 
Explicit values of I R are reported in Table |. 

Let us now define the critical crossover limit. In this case we consider the limit R — > oo, 
t — > 0, with t = R 2d /^- d H fixed. It is possible to show that 

^= ir *V(M) XH(i g ) _ /x g ) (25) 

^^^-s/M^)-,/^, (26) 

where the functions f x (t) and f^(t) are universal apart from an overall rescaling of t and 
a constant factor. Equations ( [25D and ( p6|) were predicted in Refs. [0,11] by means of a 



scaling argument and were proved to all orders in an expansion around t = oo in Ref. [^5 

The crossover functions have a well-defined behavior in the limiting cases t — and 
F-> oo. For f-> 0, Eqs. (]T7|) and (TTB|) imply 

/ x (9«j4^(1 + B"^ + ...), (27) 



x ^ 1 x 

/ € (t) « A^°r 2i/ (1 + £~t A +...). (28) 



In the limit i — > oo, for generic values of d, the crossover functions behave as 

fx® + «x^ Amf + 0(r 2A -)) , (29) 

ft(t) ~ y (1 + a^~ Amf + 0(r 2Amf )) , (30) 

where A m f = (4 — d)/2. It is important to notice that this expansion is corrected by 
logarithms whenever d — 4 — 2/k, k integer, and therefore in the interesting cases d = 2, 3. 
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For d = 3 the neglected terms in Eqs. ( p9|) and ([3(]) are of order 0(t 1 logi) and not simply 
of order O^" 1 ). A detailed derivation of these expansions and of the expressions ( p9|) and 



(p0|) is given in Ref. [35] for a much more general model than the one considered here. The 
constants a x , a^, a x , and are given, for 2 < d < 4, by 



a c = 1, 

a? = -(47r)- d/2 r(l - d/2). 



(31) 



Additional terms can be computed exactly by using the field-theoretical results of Refs. 
36TJ3Tfl , the perturbative series of Refs. [p3],|57[ , and the mean-field results of Ref. |25j, see 



Sec. IV. 



It is also possible to compute the corrections to Eqs. (j35j) and (p6|). On the basis of 
a two-loop calculation, Ref. |25j conjectured that, if the range is expressed in terms of the 
variable R defined in Eq. (||),Q then the corrections [] scale as R~ d . Explicitly, in the critical 
crossover limit we expect 



XR - f x (t) + j^h x (t) + 



(32) 
(33) 



For t — > and t oo the functions h x (t) and h^(t) have an asymptotic behavior that is 
analogous to that of the universal crossover functions f x (t) and fe(t). In Ref. [^5J the leading 
term for t — >• oo was computed, obtaining 



h x (t) w - 



(34) 



where ^ is a domain-dependent constant (see Ref. |25j for its definition). For d = 3 and 
for the domain (|3|), we have E3 w 0.058545. 



B. The fixed-length ensemble 



Given the previous results it is now a completely standard procedure [36] to obtain the 



behavior of c nj # and E\ R . For n — > 00 at i? fixed, we obtain the standard behavior 



2 This behavior can be observed only for quantities that are defined using R as scale. If we were 
considering, for instance, XR{t)p~ 2d ^ i ~ d ^ > we would of course obtain the same universal limiting 
curve, but now with corrections of order 1/p. 

3 This behavior is correct for 2 < d < 4. In two dimensions, the leading correction cannot be 
computed in a mean-field expansion and thus we do not have a theoretical prediction. However, 
numerical and theoretical arguments [18,1£ ,25] indicate that the corrections should be of order 
1/R~ 2 times powers of logi?. 
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c n ,R « C X (R) + D x (R)n~ A + •••), (35) 

El R « C E (i?) n 2 ^(l + D E (R)n- A + •••), (36) 



where 



C x (i2) 
C B (i?) 
D x (i2) 



jx(g) 

r(7) ' 

2dA g (jg)r( 7 ) 
r( 7 + 2z/) 

r( 7 - A) ' 



= r( 7 + 2z,-A) r( 7 - A) • (37) 

For R -f cx), using Eq. (P), we obtain C x (i?) -> C^R 2 ^ 1 ^^, where = A^/r( 7 ), 
with corrections of relative order R~ d . Similar relations hold for the other amplitudes. 

The critical crossover limit is trivially defined by remembering that n is the dual variable 
(in the sense of Laplace transforms) of t. Therefore, we should study the limit n — > oo, 
R — > oo with n = nR~ 2d ^ 4 ~ d ^ fixed. From Eqs. (^) and (p6) we obtain that the following 
limits exist: 

Cn,R = Cn,RPc,R ~> (38) 

<^" 8/M -^(n), (39) 

where the functions g c (n) and gE(n) are related by a Laplace transform to / x (t) and f^(t). 
Explicitly 

POO 

f x (t) = / dug c {u)e' u \ (40) 



1 f°° 

f((t)fx(t) = Yd J, du 9c(u)g E (u)e- ut . (41) 

Notice that, while the knowledge of (3 c> r is not required for the definition of ^(n), the 
critical point is needed to compute g c (n). 

The standard critical behavior is obtained for n — > oo. In this limit we have 

g c (n) « C^-\l + D™n~ A + ••■), (42) 
« C^n 2 "(l + -Dg + ••■)• (43) 

The mean-field limit corresponds to n — > 0. Using Eqs. and (3D), we obtain 

g c (n) « 1 + C c n Amf + ■ ■ • (44) 
g E {n)^2dn (l + C E n Am£ + ••■), (45) 

with corrections of order n 2Ami . In two and three dimensions additional logarithms appear. 
For d = 3 the neglected corrections to g c (n) in Eq. (HI) are of order nlogn. However, it 
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can be shown by using the field-theoretical results of App. |A 3| that the logarithmic terms 
exponentiate and that one can write 



9cin) = e^ l °z~ n g c , log (n), (46) 

where (jog is a constant and g c ,iog(n) is a function with a regular expansion in powers of n Amf 
without logarithms. The behavior of g E {n) is simpler: It has a regular expansion in powers 
of n Amf in all dimensions without logarithms ||37|| . The constants ( c and Ce can be easily 
related to a x and defined in Eqs. ( p9|) and (|3Ci|): 



Cc r(i + A mf )' (47) 

c = a x + a t a x UR] 

^ r(2 + A mf ) r(i + A mf )' 1 ] 



Using the explicit results (pTF), we obtain in three dimensions 



Cc ~ 2tt3/ 2 ' Ce ~ 6tt 3 / 2 ' (49) 

We wish now to compute the corrections to the universal crossover functions. For R — > oo, 
in the variable- length ensemble, the corrections are 0(R~ d ), see Eqs. (|32|) and (|33|), for 
2 < d < A. Thus, we expect that the universal crossover functions in the fixed-length 
ensemble have the same behavior. Therefore, we write 

c n ,R -> g c {n) + -jpk c (n), (50) 
E 2 n , R ^9E(n) + — d k E (n). (51) 
It is easy to verify by using the Euler-Mac Laurin formula that 

POO 

h x (t) = I duk c {u)e- u \ (52) 



f\(t)h x (t) + f x (t)ht(t) = — du [g E {u)k c {u) + g c (u)k E (u)} e~ ut . (53) 

The asymptotic behavior of the correction functions k c (n) and k E (ji) for n — > and n — > oo 
is analogous to that of <7 c (n) and gE^J 1 ). For n — > 0, by using Eqs. Q34]) , (|29|) , and (|30|) , we 
obtain 

fc c (0) = fe(n) = -2dE d n + O (n 1+Amf ) . (54) 



IV. FIELD-THEORY RESULTS IN THREE DIMENSIONS 

We wish now to compute the crossover functions by using field-theory methods. Consider 
the continuum </> 4 theory 
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H 



d 3 x 



(55) 



where is an N- dimensional vector — in our case N = — , and introduce the Ginzburg 
number G = u 2 ^ A ~ d "> and t = r — r c where r c is the critical value of r. Then, consider the 
limit u — ► 0, t — > 0, with tg R =t/G — tu~ 2 '^ A ~'^ fixed. In this limit we have 



X = XG - 



(56) 
(57) 



The functions F x (t SR ) and -F^(t SR ) can be computed by resumming appropriately the per- 
turbative series. There are essentially two different perturbative series one can consider: 
(a) the fixed- dimension expansion |40|j26|j27]| , which is at present the most precise one since 
seven- loop series are available p3|j57| ; (b) the so-called dimensional regularization without e- 
expansion pl|p9| , p0| that uses five- loop e-expansion results |>4]>1H . In these two schemes the 
crossover functions are expressed in terms of various renormalization-group quantities. The 
explicit expressions are reported in App. |A 1| and [A2| . For our purposes the relevant result 
is that the functions F x (t SR ) and i^(tsR) are related by simple rescalings to the crossover 
functions we have defined before f25fl. More precisely, 



/x(*) = Vx F x( st ), 



f 6 (t)= ^F^(st), 



(58) 



for appropriate constants /j x , /x^, and s. These relations are shown rigorously to all orders 



in the expansion around the mean-field limit in Ref. [^] and provide the link between 
medium-range crossover functions and field-theoretical expressions. The constants can be 
easily computed by comparing the behavior for t — > oo. In three dimensions, the functions 
F x (tsn) and i^(tsR) behave as (see App. |A 1 aj ) 



F x {t) = i[l + ^-r^ + 0{t-Ho g t)^ 



(59) 



By comparing these expansions with Eqs. 



8 = = 



and (P0|), we obtain 
1 

~ 9' 



(60) 



We can now use the explicit results of App. K| to obtain predictions for the constants 
Af>, and fl£° defined in Eqs. © and ©. We obtain in the fixed- dimension expansion 
(see App. |AT5|) 



goo 



H x XoS -7 = 0.5959 ±0.0041, 
^^s' 2v = 0.5238 ± 0.0024, 
Xis A = 2.18 ±0.18, 
£is A = 2.92 ±0.27. 



(61) 
(62) 
(63) 
(64) 



11 



For the (universal) ratio B^/B^ we obtain the more precise result 

B°° 

= 0.745 ± 0.034. (65) 

B T 

Consistent, although less precise, results can be obtained in the framework of dimensional 
regularization without e-expansion, see App. [A 2 b . 



In a completely analogous way we can derive from field theory the crossover functions 
g c (n) and gE(n). Indeed, we introduce functions G c (rit SR ) and G_E(ri SR ) in the following way: 



OO 

Ut 



F x (t) = / duG c {u)e~ u \ (66) 
Jo 

i r°° 

Fz(t)F x (t) = — J duG c (u)G E (u)e- ut . (67) 

The functions G c (^sr) an d Ge^r) can be computed perturbatively by using the cor- 
responding perturbative expressions for F x (t SR ) and Fg(t SR ). The relevant formulae are 
reported in App. [A 3| . 

In the fixed- dimension expansion, using the seven-loop results of Ref. [J57[| , we obtain 
(the six- loop result for Ge {msr) already appears in Ref. [|38|j ) 



./ 



G c (n SR ) = e 4nz log( - Kz) [1 + Az + 2irj E z 2 - 60.7295z 3 - 96.6721; 

-144.4312 5 + 2491.95/ - 5070.31/ + 0(/)] , (68) 
4 /28tt 16 n 
3 Z+ V^7~"Y 



GEinsR) = Gn : 



SR 



l + ^z+\—--^-)z 2 + 6.29688/ - 25.0573/+ 



116.135z 5 - 594.717/ + 3273.16/ + 0(/)] , (69) 
where K is a nonperturbative constant and 



1 /~ \ 1/2 
1 / ^SR X 



247T V 7T 

Explicitly 



(70) 



1 34 

\ogK = U4n 2 D 3 + -log(167r) - —, (71) 



where D 3 is a nonperturbative constant reported in App. [Ala . Numerically, using the 
estimate of D3 reported in App. |A 1 a| , we have K = 5.44(5). 



These perturbative expressions can be resummed by using the fact that the series are 
Borel summable p6|-|69| . The technical details are reported in App. |A 3| . The resummation is 
very precise for z < 1, with errors smaller than 0.2%. For larger values of z the resummation 
errors increase and the numerical integration becomes unstable: In practice we have not 
been able to compute numerically the crossover functions using Eqs. (|A87|) , (|A89|) for 
z > 5. However, in this region the crossover functions are already well approximated by 
the asymptotic expansions ( |A93| ) , (|A94| ) . The resummed expressions are well fitted by the 
following simple formulae: 
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G c (n SR ) = (1 + 50.79365z + 508.5428^ 2 + 5929.475z 3 + 10937.03,2 
G E (n SR ) = 6n SR (l + 7.6118z + 12.05135z 2 ) ' 175166 . 



4\0. 07875 



(72) 
(73) 



The expression for G^ng^) was proposed in Ref. [51] and it was obtained from a detailed 
Monte Carlo study of the Domb- Joyce model.[] We find that the perturbative results are 
very well described by these expressions, with discrepancies of less than 0.3% for z < 2. For 
larger values of z differences are slightly larger, of the order of 1%, which is, in any case, 
of the same order of the error of our resummed results. Note that the expressions (|72| ) and 
([731) exactly reproduce the small- z behaviors (j68|) and (|69| ) up to terms of order 0(z 2 ). 

The relation between the field-theory functions and g c {n), and gE{n) is straightforward. 
From Eqs. ([ 



and (|60"D we have 



g c (n) = \ c G c (pn), 



g E (n) = X E G c (pn), 



(74) 



with 



Using the results of App. 
and defined in Eqs. 



9 



A, 



9. 



(75) 



3] we can easily derive estimates for the constants C£°, 
In the fixed-dimension expansion we have 



C 



oo 7~)00 



C°° = 0.640 ± 0.005, 

A 



2.457 ±0.011, 



£>°° = 1.45 ±0.10, 

A 



5.03 ±0.48. 



(76) 
(77) 
(78) 
(79) 



If we consider the universal ratio _D°° / _D|? we obtain the more precise result 



D°° 

x 

rjoo 

E 



0.288 ±0.016. 



(80) 



We mention that from the very precise Monte Carlo results of Ref. [[HJ we would obtain 

5.57, in reasonable agreement with our results. 



2.450 and ^ 



V. ALGORITHMS 



The SAW with nearest-neighbor jumps can be very efficiently simulated by means of 
nonlocal algorithms [f7(],[7l|]. None of them can be generalized to the case at hand, and 
thus we have resorted to the dimerization algorithm (DA) |72|,73| . Although the CPU time 
needed to generate a walk increases more than any power of its length the prefactors 



4 Other representations |92-9(:.39] for Ge(tisr) are reviewed in Ref. 
stands for G B (n S R)/(6n SR ). 



]. Note that in Ref. |§ a 
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are so small that we can reach quite large lengths. It should be noticed that other algorithms 
could have probably performed better. For instance, we could have used the pruned-enriched 
Rosenbluth method of Ref. |74[ . 



Before defining the DA, let us introduce the simple-sampling algorithm (SSA). The SSA 
is the simplest algorithm for the generation of SAWs. It builds a walk recursively. Once an 
n-step i?-SAW {ujq . . . , u n } is generated, an (n + l)-step i?-SAW is obtained by choosing at 
random a new point u n+ i in D n(uj n )\{uj n } . If the new walk is self-avoiding it is kept, other- 
wise the n-step i?-SAW is discarded and the procedure starts again from scratch generating 
a new n-step -R-SAW. Since adding one step and checking for self-avoidance requires^ 0(1) 
operations, the CPU time needed to generate an n-step walk is 

n— n n 

T C PuH^^^c m /3 c m mf , (81) 

m=l 

where /3 c ,mf = 1/(Vr — !)• in the limit n — ► oo with R fixed and large, using Eqs. (|35| ) and 



(P2|), we obtain 



Tcpu(^) ~ R dj n^ +1 e anR ~ , (82) 

where a is defined by P c ,mf/Pc,R ~ 1 — aR~ d . For our model a ~ 0.035. The computer time 
increases exponentially with n although the factor in the exponential goes to zero as R~ d . 

The SSA is quite efficient in generating short walks. However, far from the Gaussian 
region it becomes too slow, because of the exponentially increasing time needed to generate a 
walk. A better algorithm is the DA | 72|j73|| . Numerically, we find DA to perform better than 



SSA for n > Vr. The DA is again a recursive algorithm. To generate an n-step walk one 
generates two n/2-step walks and concatenates them. If the resulting walk is self- avoiding, 
it is kept, otherwise the two n/2-step walks are discarded and the procedure is repeated 
again. The algorithm is recursive: in order to generate the walks of length n/2, the DA 
is used again until n/2 < n c . If n/2 < n c , we generated the walks using the SSA. In our 
implementation we chose n c m Vr. The behavior of the DA in the limit n — > oo at R fixed 



was studied in Ref. |f70| . By using the results of Sec. |HI| , one finds 

'(7-1) 



T C pv(n) ~ R qi n q2 exp 



2 log 2 



log 2 (nR 2d ^l^) 



(83) 



where q\ and q2 are exponents that depend on the specific model and on the implementation 
of the algorithm. 

Let us now discuss how to estimate E\ R and c n ^ from the simulation. Estimating E\ R 
is completely straightforward. To estimate c n> R we have used the acceptance fraction for the 
elementary moves of the two algorithms. Indeed, given an n-step i?-SAW, the probability 
of obtaining an (n + l)-step walk using the SSA is simply c n+ i[3 c ^/ c n . Thus, if we know in 



5 The self-avoidance check can be performed in a CPU-time of order one by using a hashing 
technique, see, e.g., Refs. [70,97|. 
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a given SSA simulation the number N n of generated walks of length n,we can estimate c. 
using the recursion relation 

C-n — c n-lP cm f AT , 

with the initial condition c\ = P^mS- 

Analogously, given two i?-SAWs of length n, the probability that their concatenation is 
an i?-SAW is simply C2„/c^. Therefore, if we know in a given DA simulation the number N n 
of generated walks of length n, we can compute c n using 

2 N n 

°n — C n/2 a 7" 

A n / 2 

for n > n c and then Eq. fl8~4]). 

Note that in a dimerization simulation in which we generate walks of maximal length 
^max, we obtain at the same time estimates of the observables also for a set of smaller values 
of n, i.e. for n = n max /2, n max /4, . . . These results are of course correlated, especially in 
the mean-field region where the rejection rate at each step is small. However, for our global 
observables, the correlation should be small. Analogously, when we use the SSA, we can 
compute the observables for all values of n, although the results are in this case strongly 
correlated. 



(84) 



(85) 



VI. NUMERICAL RESULTS 

We have performed an extensive simulation using three-dimensional walks with n < 
66560 and 2 < p < 12. Notice that the values of p are particularly large: for p = 12, in 
the spin language, each spin interacts with 2624 neighbors. The advantage of working with 
SAWs is the absence of finite-size effects — we work in the infinite-volume limit — and the 
possibility of reaching large values of the correlation length.^ The raw data for the largest 
valuesQ of n and several values of R are reported in Tables |II| and 



6 In order to compare with spin models, it is interesting to relate our values of R with £r(/3), 
where, for each n, we consider the value of (3 such that (n)p = n. In the critical limit we find the 
relation 



^> ~ 6 r( 7 ) 7 2- E »,* 



Therefore, we reach f » 35,80,300,500,350,500,500,700,700,500, respectively for R = 
2,3,4,5,6,7,8,9,10,12. 

7 For p = 2, 7, 12, we have computed the observables for additional values of n, larger than the 
minimum value of n reported in the Table. More precisely, by using the SSA, we have computed 
the observables for all n < 30, 400, 3000 respectively. 
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In Sec. [VI A| we will determine /3 Cj r from our numerical data and we will explicitly check 
the theoretical predictions for the large- R behavior of /3 C) # of Ref. [25j presented in Sec. 
1 11 A| . In Sec. |V1 B| we will compute the critical crossover functions and we will compare 
them with the field-theoretical results of Sec. TV . 



A. Determination of (3 Ci r 

In order to compute (3 Ci r we define 

Cn ,R 



(3ett,R( n ) 



_9c,th( n ). 



-l/n 



(86) 



where g c ,th{n) is the theoretical crossover function: in our numerical determination of 
PeB,R(n) we will use0 Eqs. (|74]) and (|A106J) . By using Eqs. (|35|) and fl42"D , we obtain 
for n — > oo 



(3eK,R{ n ) = A 



c,R 



n 



log 



C X (R)R 6 ^ 



l+A 



n 



(D X (R) - D°°R 6A ) + ■ ■ 



(87) 



Using the asymptotic expansions ( pT|) and the relations (p7|), we obtain for i? — > oo 

?6A-3 



1 + _^!L _l a2R 



nR 3 



n 



l+A 



+ 



(88) 



where a% and «2 are /^-independent constants. This expression shows the advantage of the 
definition ( p6| ) over the common one in which one simply considers (c n ) _1//n . Indeed, with 
our choice, the l/n correction vanishes for R — >• oo while the l/n 1+A remains approximately 
constant (A rs 1/2); with the other one, we would have corrections of order log(nR~ 6 )/n 
and R 6A /n 1+A . This improved behavior is particularly important, since for large R we are 
quite far from the Wilson-Fisher point, and thus a reduction of the scaling corrections is 
essential in order to obtain precise estimates of (3 C ,R- I n order to determine /3 C) _r we have 
performed fits of the formf] 



/3eS,R( n ) V R = Pc,rVr + - + ~^A- 



n 



(89) 



8 In principle, we could have used Eq. ([72]) for G c (z). However, here it is important to use an 
expression for <? Cj th(^) which is precise in the Wilson-Fisher region. For this reason, we have chosen 
to use Eq. §KW§ . 



9 In the fits we also used estimates obtained by means of the SSA. In order to reduce the correlations 
and to avoid the fit to be dominated by data with small values of n, we have used only some of 
them. We considered the values corresponding to n = no [ e ( fc + fc o)/5j ; ^ g JV", where [^J is the 
largest integer smaller than x and no and ko are integers. 
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assuming A = 1/2. We have repeated the fit several times, considering each time only the 
data satisfying n > n min . The final results, reported in Table J^, correspond to the smallest 
n min for which x 2 /d.o.f. ~ 1 (d.o.f. is the number of degrees of freedom). Notice that, by 
rescaling b by R 3 , the coefficients a and b should become ^-independent as R increases. This 
is evident for a and indeed we can roughly estimate a ~ —0.015(5) for R — > oo. This allows 
us to compute the leading correction to C X (R). We have 

C X (R) « C~i? 6{1 - 7) (1 + k x R~ 3 + ■■•), (90) 

where The results for b are less stable, but still reasonably compatible with a 

constant for large R. In order to understand the systematic errors due to the truncation 
(|S9|) we have repeated the fit with an additional correction: 

fcz, R {n)V R = (3 C , R V R + - + ^ + SfL. (91) 

The results for (3 C ^ do not differ significantly from those of the fit with c = 0, except for 
p — 12, where the difference is approximately three combined error bars. Therefore, our 
final estimates should be quite reliable. As an additional check we have compared our results 
with the theoretical prediction (1221) . If we indicate with /^jj the expansion (22) neglecting 
terms of order R~ 8 , we define 

- # i% - 1 ) • (92) 



A 



c,R 



If we have correctly determined (3 C! r, r 2j efr should converge to the constant r 2 as R — ► oo, 
with corrections of order log R 2 /R. The plot of t 2)C r is reported in Fig. [l|. A fit of the form 
i~2,efr = r 2 + b/R gives 

r 2 = -0.00814(21), (93) 
b = 0.00277(19), including all data with p > 3. 



B. Determination of the critical crossover functions 

We wish now to determine the critical crossover functions in three dimensions. We begin 
by studying the function c ni R, cf. Eq. fl3~8"l). The function is reported in Fig. |2] (upper 
graph) together with the theoretical prediction obtained by using Eqs. ( ff4[) and (|72|) . Note 



that there is no free parameter in the theoretical curve. We observe a very good agreement, 
especially in the Wilson-Fisher region. Systematic deviations are observed for smaller values 
of n. In order to understand the role of the deviations we report in Fig. (lower graph) the 
same data, but now we exclude all points with n < Vr/2. The agreement is now perfect for 
all p > 3. We thus clearly see that the crossover behavior requires n>i? 3 . In particular, 
mean-field behavior is always observed for n -C -R 6 if R is sufficiently large, but it is described 
by the critical crossover curves only if R 3 n R 6 . 
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In order to see the corrections to scaling, in Fig. we report the ratio c n ,n/ g Ct th{n) which 
should converge to 1 as R — > oo. Corrections to scaling are clearly evident, points with 
different values of R lying on different curves that indeed converge to 1 as R — > oo. These 
corrections are predicted to scale as R~ d . To check this behavior we considered 

A C ;n,H = R 3 f-^r " l) , (94) 
\9cM n ) J 

which should converge to k c (n)/g c (n) in the crossover limit. Using the expected asymptotic 
behavior of k c (n) and g c {n), A c;n> # converges to a constant both for n — > and n — > oo. For 
n -► 0, k c (p)/g c (0) = -E 3 « -0.059. 

The numerical results are reported in Fig. |], where the error on A c . n r has been computed 
by considering the error on /3 c> r, c n ^, and on the theoretical curve.0 A reasonably good 



scaling is observed, confirming the results of Ref. [2l|. Also the prediction k c (0)/g c (0) 
—0.059 is fully compatible with the data. 

In order to perform a more precise check, we have also considered the quantity 



Q n . R = (95) 

C2n.R 



that converges to g c (n) 2 / g c (2n) in the critical crossover limit. For Q U) r we do not need the 
value of (3 c> r and thus a source of error is avoided. In Fig. |5] we show Q n ,R together with 
the theoretical prediction and 

A Q;n ,* = i? 3 ( 9^9c^ _ A (96) 
V 9cM n ) J 

The errors on Aq^^r have been computed as we did for A c;ni #. The agreement between 
the numerical data and the theoretical prediction is very good. Also, A.Q ]nj R shows a nice 
scaling behavior confirming that the corrections scale as R~ d . 

Let us finally discuss the effective exponent j e ff(n,R). A standard definition would be 

7 eff (n, J R) = 1 + n . (97) 

an 

However, this definition is not easy to use in numerical simulations since it involves the 
derivative with respect to n. Here, we will use the definition 

leS (n,R)= 1 + J-log(^), (98) 



log 2 \ C nt R 

which interpolates between the SAW value 7 = 1.1575 and the mean-field value 7 = 1. The 
results are reported in Fig. ^| together with the theoretical prediction. The agreement is 



10 We assume g c ,th( n ) to have an error of 0.3% for n > 1, and of 0.3n % for n < 1. The smaller 
error for n < 1 is due to the fact that our representation (72) has the exact leading behavior for 
z -> 0, cf. Sec. HV. 
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very good. Note that in the Wilson-Fisher region the numerical data are well approximated 
by the field-theoretical prediction only for R > 4. For smaller values of R, the corrections 
are important, as it was already noticed in the Ising model simulations • 



The analysis we have performed for c n can be repeated for E\ R . In this case however 
the errors are smaller since the critical crossover functions do not depend on f3 c ,R- In Fig. [7] 
(upper graph) we report our results for E% R together with the prediction gE ,th(n) obtained 
by using the field-theory result ( ff3|) and the relations ([74]) , (|75|). Note that there is no free 
parameter in the theoretical curve. The agreement is very good although one can see clearly 
the presence of corrections to scaling. 

We wish now to compute the correction curve fc^(n). For this purpose we consider 

A w = R 3 I -%r - l) , (99) 

\9EM n ) J 

that converges to &E;(n) / 9E(n) as R — > oo. The plot of Ae-u,r is reported in Fig. |7| (lower 
graph), where we have taken only into account the error on E\ R . A good scaling behavior 
is observed confirming the theoretical prediction for the corrections. Moreover, this nice 
scaling behavior is also an indication that the approximation (|73|) can be considered at our 
level of precision practically exact. Note also that the prediction fcE(0)/<7e(0) = —0.059, cf. 
Eq. ( f>4|) , is in good agreement with our data. 

Let us finally discuss the effective exponent v e $. As we did for 7 eff , instead of 

p eS {n,R)=n ° n ' R , (100) 
an 



we will consider 



/? 2 

Cj 2n,R 



^<"'*>= 21^2 ^ iff 1 ' <101) 

that is easier to compute numerically. A graph of this quantity is reported in Fig. || It 
shows the expected behavior: for n — > it converges to 1/2, while for n —>■ oo it converges to 
v saw ~ 0.588. The agreement with the perturbative prediction is quite good in the random- 
walk region. On the other hand, as n increases the corrections increase, in agreeement with 



similar results obtained for the Ising model pO|pl 



VII. CONCLUSIONS 

In this paper we have studied the critical crossover limit for a model of walks and we have 
verified numerically the following statements, predicted by field-theoretical and mean-field 
methods 

1) The critical crossover functions in medium-range models coincide with the field- 
theoretical crossover curves. The nonuniversal constants can be determined by com- 
puting the corrections to the mean-field behavior. 
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2) The asymptotic behavior of (3 Ct R for R — > oo can be computed exactly in lattice models 
up to corrections of relative order R~ 8 by determining the first corrections to the mean- 
field limit and exploiting the field-theoretical model. 

3) The corrections to the critical crossover functions decrease as R~ d once R is defined 
as in Eq. <^). 

Our numerical results can also be used to determine phenomenological expressions that 
describe the data for all values of R. Here, different procedures can be used. One can 
consider the phenomenological model of Ref. |75j (with the modifications discussed in Ref. 



25| to make it compatible with the theoretical predictions) or use the procedure proposed 



in Ref. f2~5fl . The idea is to write 



c n ,R = 9c{n) + R 3 k c (n), 
El R = g E (n) + R- 3 k E (n) } (102) 

and use a simple parametrization for the correction terms. Here, we approximate 

-0.059 + an 1 / 2 + bn 
k[n) = g{n) — - — - — 103 
1 + cn L i A + an 



both for k c and k E , where we have used the asymptotic behavior fl5J). The parameters a, 
b, c, and d are determined by fitting the numerical data. The best results are obtained for 
a = -61, b = -1.06, c = 1830, d = 87 (function k c ) and a = -23, b = 0.8505, c = 972, 
d = 32 (function k E ). These fitting functions provide phenomenological expressions that 
correctly describe our data for all values of p. The corresponding effective exponents are 
reported in Fig. |9| and show the typical behavior that has been found in simulations of the 
Ising model. We have also included in the figure the curves corresponding to p = 3/2, to 
show that for p small the phenomenological expressions show a nonmonotonic behavior that 
is not present in the critical crossover curve. 

Although our main motivation was the understanding of spin models, the results of 



this paper are also relevant in the context of polymers. Indeed, as is well known [|35| , |36 
noninteracting SAWs describe the universal behavior of homopolymers in dilute solutions 
above the O-temperature. In the polymer context, however, it is more interesting to consider 
a different model with medium-range interactions. Supposing, for simplicity, to be in the 
continuum (off-lattice), we can define a SAW in the following way. A SAW with medium 
range interactions is a collection of point {u , • • • u>i G R d , such that \u>i — oJ i+ i\ = p 

and \uji — ujj\ > a for all i ^ j. In this case the relevant scale is p/a and the crossover 
limit is obtained for p/a — > oo, n — > oo, with n = n(p/a)~ 2d ^ 4: ~ d ^ fixed. For this model, the 
critical crossover functions can also be computed using Eq. (|7"4f ), although with different 
nonuniversal p-independent constants X E , A c , and p. Thus, the results presented here are 
relevant for the description of polymeric systems in which the macromolecular persistence is 
much larger than the molecular scale. In practice, we expect the description to be reasonably 
accurate when p/a > 3. 
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APPENDIX A: CRITICAL CROSSOVER FUNCTIONS FROM FIELD THEORY 



In this Appendix we will compute the critical crossover functions for the polymer case 
using field-theory methods. In Section |A 1| we will use the approach of Refs. p6| , p7| , while 
in Section |A 2j we will present the results obtained using the method of Refs. l|4l,29,[3~0|. The 
first approach provides the most precise estimates and it will be applied in Sect. |A 3| to 
obtain numerical results for polymers generalizing Ref. [pq| . 



1. Crossover functions in the fixed-dimension expansion 



a. General results 



In this Section we report the critical crossover functions using the approach of Refs. 
26| , |27[| . We start from the expressions for F x (t) and F^(t): 



F x( f ) = X* exp 



? , dX v(x)W(x)_ 



Ft® = (t) 2 exp 



dx 



1 



W(x) 



(Al) 
(A2) 



where t is related to the zero-momentum four-point renormalized coupling g by 



dx 



7(ar) 



v(x)W(x\ 



exp 



dz — — — 

, »{z)W{z) 



(A3) 



7(x), z/(x), and W(x) are the standard renormalization-group (RG) functions, g* is the 
critical value of g defined by W(g*) = 0, and x*i £*> ^o, and yo are normalization constants. 

The expressions (|Al]), (|A2|) , and ( |A3|) are valid for any dimension d < 4. The first two 
equations are always well defined, while Eq. ( |A3[ ) has been obtained with the additional 
hypothesis that the integral over x is well defined when the integration is extended up to 
g* . This hypothesis is verified when the system becomes critical at a finite value of /3 and 
shows a standard critical behavior. In our case, = 0, this is true for all 1 < d < 4. 

We normalize the coupling g as in Refs. |63| , |76| so that in the perturbative limit g — > 0, 
t — > oo, we have 



g « — i_ r (2--) = A d ? d - 4 )/ 2 

y 3(4vr) d / 2 V 2 J 

This implies that for y — > we have t m (yo/Xd) 2 ^ d ~^ and (C) 2 t 
normalization, the previous equations can be written as 



X*t 



F x {t) = (6th?) 2 exp 



Fz(t) = (Girg) 2 exp 



dx 



-2 / dx 
'0 



7(x) 



+ 



v(x)W(x) x 



1 1 

W(x) x 



(A4) 
1. With this 

(A5) 
(A6) 
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and 



1 



(67rs0 2 

i r 9 dx 



(67r) 2 J x 2 [u(x)W(x) 



7(x) 



exp 



dz 



v{z)W{z) z 



where D 3 is a nonperturbative constant given by 



1 



1 



(6tt) 2 
1 f 9 * dx 



3 1, * 



(g*) 2 2g* 4 
7(x) 

exp 



(67r) 2 7 x 2 [z/(a;)W / (x) 



dz 



LJO 



v(z)W(z) + 2 



2 3 x 
x ~~ 2 ~ 4 



+ 



4 



(A7) 



(A8) 



Numerically, D3 = 0.002473(6). Expressions for a general iV- vector model and in two 
dimensions can be found in Ref. . 

For t — > 00, in three dimensions we obtain 



F x (t) = ? 



- — r 1 / 2 L_i og ( 367r 2 t) 

12tt 288vr 2 t V 7 

59 1 D 3 



2592tt 2 t t 



^+^ + 0(r 3/2 logt) 



(A9) 



1 + 



1 



127T 
11 



-t 



-1/2 



1 



1 £>3 

'486tt 2 T + T 



dog (367T 2 t) 

2887r 2 t V ; 
+ 0(r 3/2 log?) 



(A10) 



6. Asymptotic behavior near the Wilson-Fisher point 



Let us now compute the asymptotic behavior of the crossover functions for t — > 0. This 
requires the determination of the expansion of the various RG functions in the limit g — > g*. 
As it has been extensively discussed in the literature [ ^0| , [r^ , [79| , [7^ , ^7| , ^0| -p4| , these functions 
are singular at g = g*. General arguments predict a behavior of the form |79|J8"2| 



W(g) = -u{g* -g)+ Wl (g* - gf + w 2 (g* - g)^ A + w 3 (g* 

1(9) = 1 + 11(9* -9) + 12(9* - gf + 1,(9* - g) A2/A + 74(0* 



9) 



v(g) = " + vi(g* - 9) + Mg* - gf + Mg* - g) 



A 2 /A 



Mg*-g) 1/A + 



1+1/A 

,)VA 
1/A 



+ 



+ ■ ■ 



(All) 
(A12) 
(A13) 



This nonanalytic behavior makes the determination of the corrections extremely difficult. 
For instance, since one expects A 2 /A to be close to 2 |85| , it is practically impossible to 
determine Wi and w 2 in the /^-function since these two terms are essentially degenerate. The 
only subleading coefficients that can be reliably determined are 71 and v\. Indeed, since 
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A f» 1/2 and A2/A pa 2, the next-to-leading correction behave as (g* — g)~ 2 , so that ^f'(g) 
and v'{g) should be reasonably smooth for g — > g*. 

The computation is straightforward and we only report the final results. The crossover 
functions can be expanded for t — > as 



F x (t) = Xo* -7 [l + Xit A + X2t + X 3 t A2 + X 4t 2A + -- •] 

F c (t) = &~ 2v [1 + + + 6J Aa + e 4 ? A + •••]. 

We obtain for the leading term and the first correction: 



dx 



Xo = (Qrrg*) 2 ^' exp 
e 2 =(67r^) 2 P^exp|-2^ dx 



70) 



2 

+ - + 



7 



v{x)W(x) x A(g* — x) 

1 1 1 
+ - + 



W(x) x oo(g* — x) 



Xi 

6 



9 1i . f _ A 



A(l + A) 



t -A [71/1(1 + A) - 1/71A] 



7A(1 + A) 



(A14) 
(A15) 



(A16) 

(A17) 
(A18) 
(A19) 



where 



7 



(Qirg*) 



exp 



1 2 

u(x)W(x) X A((yf* 



(A20) 



In order to estimate these constants we have used the seven-loop results of Ref. |57 and the 



resummation technique of Ref. ||76|| . Errors due to the resummation have been determined by 
using the general procedure of Ref. 
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In order to compute xo an d Co we have performed 
resummations keeping g*, 7, u, and A as free parameters (the dependence on A cancels in 
Xo and £q if one uses the explicit expression for t) . We obtain finally 

Xo = 0.4216 ± 0.0006 - 4( 7 - 1.1575) - O.lfo* - 1.395), (A21) 
£0 = 0.3565 ± 0.0002 + 0.4(7 - 1.1575) - 0.1(p* - 1.395) - 8(1/ - 0.58758), (A22) 
t= [0.99 ± 0.004 + 0.9(7 - 1-1575) - 8(#* - 1.395) - 9(A - 0.515)] • 10~ 3 , (A23) 

where the first error is related to the uncertainty in the resummation, while the other terms 
indicate the variation of the estimate with changes in the values of the critical exponents and 
of g*. To obtain the final results, we must decide which estimates to use for 7, u, A, and g*. 
In principle, we could use the values that have been determined from the resummation of the 
perturbative expansion in three dimensions [E3], i.e. g* = 1.413±0.006, 7 = 1.1596±0.0020, 
v = 0.5882 ±0.0011, and A = 0.478 ±0.010. However, we believe the Monte Carlo estimates 
of the critical exponents to be more reliable, and thus we have used the values reported 
in Eq. (p0|). The field-theoretic estimate of g* is probably also not reliable since it differs 
from the estimates obtained using different methods. Indeed, from the e-expansion one 
estimates g* = 1.396 ± 0.020 ||86|| , while the extrapolation of exact-enumeration series gives 
g* = 1.388 ± 0.005 p7l: One observes a systematic discrepancy, which we believe to be 
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due to the nonanalytic structure of the /3-function that is not properly taken into account 
in the analysis.^ This problem should also appear in our analysis since we use the same 



resummation technique. In the reanalysis of Ref. of the perturbative series for the 
exponent 7 it was shown that the systematic error could be reduced, obtaining field-theory 
estimates in close agreement with the Monte Carlo results, if one uses g* ~ 1.395. Therefore, 
we use g* = 1.395±0.015, where the error is such to include all estimates. Our final estimates 
are 

Xo = 0.4216 ±0.0029, (A24) 
£2 = 0.3565 ±0.0016, (A25) 
t = (0.99 ±0.21) • 10" 3 . (A26) 

To compute %i an d £1, we have analyzed the series of the derivative of j(g) and v(g). We 
obtained 

7l = -i(g*) = -0.1071 ± 0.0013 + 0.007(#* - 1.395), (A27) 
Vl = -i/(g*) = -0.0659 ± 0.0018 - O.Ollfa* - 1.395). (A28) 

The ratio 71/1/1 has already been computed in Ref. |57j finding 

^ = 1.31 ±0.05- 1.7(0* - 1.39) (A29) 

that however differs significantly from our result 

— = 1.62 ±0.05 -0.4(0* - 1.395). (A30) 



Using the estimates of the critical exponents reported in Eq. ( p0|) and, as before, g* = 
1.395 ± 0.015, we obtain 

Xi = 6.8±0.8, (A31) 

£i = 9.1 ±1.1. (A32) 

Notice that a significant fraction of the error is due to the uncertainty on t. The error is 
largely reduced if we consider the ratio Xi/£i- We obtain 

§i = 0.745 ± 0.034. (A33) 
£1 

Notice that, if we use the estimate ( |A29|) we would obtain Xi/£i ~ 0.56. The ratio Xi/£i 
can also be computed in e-expansion using the 0(e 2 ) series of Ref. 



£1 8 V 12 256 

where A ~ 1.171854. We obtain Xi/£i = 0.85 ± 0.10, where the error is purely indicative 
because of the shortness of the series. 
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A different resummation that tries to keep track of the nonanalyticities of the /3-fimction give 



[57] g* « 1.39. 
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2. Crossover functions in dimensional regularization without e-expansion 



a. General results 



In this Section we will study the critical crossover functions using the minimal renormal- 
ization scheme proposed in Refs. |4l|j29| . |30[| . We start from the expressions [^9],[30| 



F x (t) = X *F(u)~ l exp 



1M{X) 



dx 

u M (x)W M (x) 



F&) = (C) 2 exp 



-2 



dx 



W M (x) 



(A35) 
(A36) 



where t is related to the minimal-subtraction renormalized coupling u by 

1 



t = -to 



dx ~~~— t~t exp 

W M (x) 



dz 



vm(z)W m (z) 



(A37) 



Here jm(x), Um(x), and Wm(x) are the standard RG functions computed in dimensional 
regularization, P(x) and F(x) are functions defined in Refs. p9],|30| that will be explicitly 
given below, and u* the critical value of u defined by Wm(u*) = 0. The constant %* i V, to, 
and Ui are obtained by requiring that, for t — > 0, F x (t) and F^(t) behave as in Eq. (^). The 
RG functions have been computed to five-loop order in Refs. H^JH^]. Explicitly we haveff] 



W M {u) = W M {u,l), 

2 + r] 2 (u) 
Jm(u) = [2 - 773 (¥)] u M (u), 



(A38) 
(A39) 
(A40) 



where 
W M (u,e) 



en, _L 4„.2 _ 7„,3 1 2960+2112C(3) 4 , -196648+2816tt/5-223872C(3)-357120C(5) 5 , 
eU-^ 3 U 6 1728 U 62208 " 



13177344-67584tt 4 - 3174407T 6 /21+21029376C(3)+2506752C(3) 2 +42261504C(5)+59383296C(7) 6 

1990656 U 



+ 0{u 7 ), (A41) 



3 « -r 36 u 144 u 



-31060-3527r 4 /5-3264C(3) 4 



62208 



3166528+426887r 4 /5+396807r 6 /21+1528704C(3)-446976C(3) 2 +55296C(5) 5 , /-)/ 6n 

2985984 M + U l U J' 



(A42) 



„ ( v ) - 1 ,.2 _ 1 ..3 , 125 4 , -77056-1408^/5+8832C(3) 5 , 6>> 



(A43) 



To compute the functions P(u) and F(u) let us first recall the relation between the bare 
coupling mq and the renormalized coupling u 



12-1 



We normalize u as in Ref. |p8[ . Our coupling u differs from the coupling used in Refs. [2S.3C] 



by a factor of 4! and from the definition used in Refs. [64,65| by a factor of 2. 



25 



u = uA 3 1 uZ u (u)Z ip {u) 



Zuiv^Z^iu) 2 = exp 



du' 



1 



W M (u> 



2 
u' 



(A44) 
(A45) 



Here \x is the renormalization scale and Ad is a constant that depends on the specific renor- 
malization scheme. Of course physical results should not depend on it. This fact can be 
easily verified noticing that Ad can be absorbed in the definition of fi, and that, by con- 
struction, all physical quantities are independent of u. However, different choices of Ad give 
rise to different perturbative series providing different results at the intermediate stages of 
the calculation. This freedom may be used as a further check of the real uncertainty of the 
final results. In Refs. PH^OJ the authors use 



A d = S d T(3 - d/2)V(d/2 

a 2 
d {4n) d / 2 ^{d/2) , 



(A46) 



a choice that makes the one-loop corrections vanish in many observables. However, in order 
to understand the size of the systematic errors, we will also use Ad = Sd and Ad = 45*^. The 
functions F(u) and P(u) are obtained from 



F(u) = Z^(u)F (A 3 uZ u (u)Z^(u 



Z v (u) 



exp 



du 



W M (W 



and 



P{u) = Z^(u)- l P (A^uZ u (u)Z^u)- 2 ), 

/ V2{u') 



Z^i\u 



exp 



du 



W M (u> 



(A47) 
(A48) 

(A49) 
(A50) 



where F (x) and Po(^) can be derived from the five-loop results of Refs. [p7| , p0| 



F (x) 
Po(x) 



1 + 



1 

7776tt 



-X 



8.83291 • 10~ 7 x 3 + 6.17241 • KT 8 ^ 
4.73993 • KT 9 x 5 + 0(x 6 ), 



i i J_ r i 

x ~ 24tt 288tt 2 



x 2 + 1.81785 ■ 10~ 5 a; 3 - 1.24518 • 10~ 6 x 4 



+ 1.01097- 10" 7 x 5 + 0(x 6 ). 
When Ad is given by Eq. (|A46 ) we have 
F(u) 



(A51) 
(A52) 



1 - yI^u 2 - 5.52510 10~ 3 u 3 - 4.01633 10~ 3 u 4 - 1.92954 10" 3 w 5 + 0(u 6 ), 



P(u) = 1 - \u + ±u 2 - 0.0318279 u 3 + 0.0326664 u A - 0.0644065 u 5 + 0(u G 



(A53) 



(A54) 
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b. Asymptotic behavior near the Wilson-Fisher point 



In order to compute the asymptotic behavior of the crossover functions for t — > 0, we 
need the expansion of the various RG functions in the limit u —>■ u* . It has been argued 189 



that the scheme we are presenting has an important advantage over the approach described 
in App. [A 1| : The RG functions are expected to be analytic at the critical point. The 
reason is that the RG functions are essentially dimension- independent, while u* depends on 
e being the solution of Wm(u*, e) = 0. Notice, however, that this argument does not exclude 
the presence of singular terms for u = since this scheme is essentially four-dimensional. 
Answering this question is in any case nontrivial since it requires a nonperturbative definition 
of the RG functions. In the following we will assume the following analytic expansions: 



W M (u) = -uj(u* - u) + wi(u* -u) 2 + ... 
Jm(u) = 7 + Ji(u* -u) + ... 
Vm{u) = v + v\{u* —u) + ... 
F(u) = r + f 1 (u*-u) + ... 
P{u) = p* + pi(u* — u) + . . . 



(A55) 
(A56) 
(A57) 
(A58) 
(A59) 



Expanding the crossover functions for t — > as in Eqs. ( A14 ) and (|A15 ), we obtain for the 
leading term and the first correction: 



Xo 



e 



[KU 



*\2 



-t 1 exp 



r 



[ku*) 2 ^ exp < -2 



dx 



1MK.X) 



2 

+ - + 



7 



dx 



u m (x)Wm(x) x A(u* — x) 
1 



Xi = -u 



1 2 

+ 



Wm(x) x lo(u* — x) 



fi 7i 7_ 

/* A 1 + A\ u 



"1 , Pi 
p* 



t * 

£i = u 



2u 



Pi "I \£-A 



1 + A \p* Ais 



(A60) 
(A61) 
(A62) 
(A63) 



where 



2p*v 

1 — 7\2 ex P 

[KU V 



a* 



dx 



u m (x)Wm(x) x A(u* — x) 



The normalization k is related to the choice of Ad in Eq. (|A45|) by 



K = 4vr(l + 2g 3 ), 

S d T(3 - d/2)T{d/2 - I) 



(A64) 



(A65) 



A, 



Qd 



2(4 - d)A B 



The estimate of critical quantities requires a resummation of the perturbative series. For 
this purpose we will use the large-order behavior of the coefficients |9^j9l| given by 



c k ~k\ (-|) fe A; 6 [l + 0(l/A;)] 



(A66) 
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and the numerical method of Refs. [|76| , |82 . 

We have first of all determined u*. We obtain u* = 1.1 ± 0.1. Although this result is 
consistent with the estimate of Ref. |[29|| , u* = 1.092 ± 0.012, the error bar is much larger. 
However, 



82 



, using our algorithm 
be so small. On the other hand 



, we have been unable to understand how the error can 
as we shall see later by comparing our results with the 
estimates of the previous Section and by checking their independence on Ad, our error bars 
look reasonable and at most overestimated by a factor of two. In the following we report 
various estimates keeping u* as a free variable. We have 



7i — - 


-U.14U It u.uiu — U.Ul^xi • 


- l.lj, 




v\ = - 


-0.089 ± 0.010 - 0.02(w • 


- LI), 


(A68) 




' 0.991 ±0.004 - 0.02(w* 


- LI), 




r = < 


0.987 ± 0.006 - 0.02(w* 


-l.l), 


(A69) 




k 0.983 ± 0.006 - 0.02(m* 


-l.l), 






' 0.932 ± 0.006 - 0.05(w* 


- LI), 




P * = < 


0.825 ± 0.015 -0.15(u* 


-l.l), 


(A70) 




^ 0.670 ± 0.040 - 0.30(w* 


-l.l), 






' 0.026 ± 0.007 ±0.03(w* 


-l.l), 




... 


0.039 ± 0.012 ±0.04(w* 


- l.l), 


(A7l) 




^ 0.051 ± 0.018 + 0.06(ti* 


-l.l), 






' 0.079 ±0.015 + 0.03(w* 


-l.l), 




pi = < 


0.182 ± 0.012 ±0.02(w* 


-l.l), 


(A72) 




^ 0.399 ± 0.030 + 0.13(u* 


-l.l). 





For /*, p*, fi, and p\ we report three estimates corresponding to A 3 = S3, 1/47T, 4Ss re- 
spectively. Notice that in most cases the uncertainty on u* is negligible compared to the 
resummation errors. We obtain finally 



Xo = 0.420(15), 0.422(19), 0.447(42), 
Co = 0.358(15), 0.363(19), 0.395(38), 
X i = 6.7(4.2), 6.5(4.1), 5.4(3.4), 
fi = 10.1(6.3), 9.2(5.7), 6.7(4.2), 



(A73) 
(A74) 
(A75) 
(A76) 



where the three different estimates correspond to A 3 = S3, 1/47T, 4Ss respectively. As before, 
a more precise estimate is obtained if one considers Xi/£i- We obtain 



= 0.66(13), 0.71(10), 0.81(7). 



(A77) 



As expected, these results are independent of the value of A 3 within error bars. Notice that 
the difference among the estimates of the same quantity is of the same order of the error 
bars, thereby confirming the correctness of our error estimates. The final results are also 
in good agrement with, although less precise than, the results presented in the previous 
Section. Notice that the estimate of obtained here is compatible with the estimate 

obtained in Sec. [A 1 b|. Instead, the result of Ref. p7[ is somewhat too small. 
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3. Polymer critical crossover functions 



Let us now compute the critical crossover functions in terms of n. We will only consider 
the approach described in Sect. [A 1[ since it appears to be the most precise one. From Eq. 
(HO) we have 



G c (n) 



c+ioo 7, 

— e nt F (t) 
■ 2m x{ h 



(A78) 



Changing variables from t to g we obtain 



c 



2m 



"(9)W(g) 



exp 



(A79) 



where C may be taken as a circle of the form g = R(l + e -4<? j 
satisfying < R < g*/2. 

To compare with the results of Ref. PS) we introduce their notations 



l-> hi i y{z)W{z)\ ' 

-7T < (b < 7T, with i? fixed 



7(fl0 



j(<?) = 2(1 



1 



70), 

(l- 7 )/A 



= exp I J dx 

"9 



exp 



dx 



7-1 



2-?'(x) 2 

r -^- H 1 

(3{x) x A(g* — x) 



_(3{x) A(g*-x) 
1 



E (g) 



dx ( 1 

9 * V 9* 



1/A 



A(ar) 
a; 2 /?(a;) 



where 7 and A are the standard critical exponents. Then, we obtain 

dg J(g) 



GJn) 



c 2mP{g) 



exp 



n 



18m 



E (g) 



In order to have the same definitions of Refs. |37lJ38|l , let us also introduce 



n 



24vr 3 / 2 ' 



We obtain finally 



GJz) 



J (d) e 32nz 2 E (g) 

c ™ /3{g) 



(A80) 
(A81) 

(A82) 
(A83) 

(A84) 



(A85) 



(A86) 



(A87) 



Expanding the previous expression for z — > 0, we obtain the perturbative expansion 
The computation of Ge{z) is analogous. Starting from 
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G E (n) 



GJn) 



c+ioo j j. 

— e^F x m(t), 

c-ioo 



(A88) 



we obtain for the swelling factor |38| Se(z) = Ge(z)/(6ti) 



Se(z) 



(Id J(g)K(g) 32^E (g) 



16ttz 2 G c (z) J c ttz 0(g)E(g) 



where 



Kg) ■■ 

K{g) 

E(g) 



(Mg) : 
V i(g) 



g_ 

9* 



exp 



dx 



k(x) 7 — 2v 



g' 



g_ 

g* 



t/a 



exp < 2 / dx 



j3(x) A(g*-x) 
7 



1 1 
+ - + 



(3{x) x 2A(g* 



x 



(A89) 

(A90) 
(A91) 
(A92) 



From the previous expressions we can compute the asymptotic behavior of G c (z) and Se(z) 
for z — > oo. We have 



G c (z) = g c0 z 2 ^ (1 + g cl z- 2A + g c2 z^ + g c ,z~ 2 ^ + . . .) 



SfXz) 



seqZ 



Au-2 



(l + s E1 : 



-2A 



+ S E2 Z + S£ 3 Z 



-2A 2 



+...)• 



(A93) 
(A94) 



These expansions can be related to the expansions of F x (t) and F^(t) for t — *■ 0. Using the 
results of App. |A 1 b , we have 



g c o 

seo 

ga 

SEl 



(24vr 3 / 2 ) 



Xo 

r(7) 

am 

r( 7 + 2z/ 

xir(7) 



3/2\ 2(7-1) 



2.117 ±0.020, 



(24tt 3 / 2 ) 



3/2^ 2 (2^-l) 



1.549 ±0.007, 



r( 7 - A) ' 
(xi±£i)r(7±2z/) 



-2A 



= 0.029 ±0.005 

Xi r (7) 



(24vr 3 / 2 ) 



-2A 



0.101 ±0.016. 



r( 7 ±2z/-A) r(7-A) 
If we consider the ratio g c \j s e\ the error is largely reduced and we have 

— = 0.288 ±0.016. 

SEl 



(A95) 
(A96) 
(A97) 
(A98) 

(A99) 



The constants seo and s E i have already been computed in Ref. |38[, finding seo ~ 1.5310 
and sei ~ 0.1204, in reasonable agreement with our results. The constants seo and s^i have 
also been determined by a Monte Carlo simulation of the Domb- Joyce model [51], obtaining 

(A100) 



seo = lim B R (u) « 1.54654 



s E i = lim b R (uj) 

LO— >0 



0.11498 



(A101) 
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where Br(lj) and 6r(u;) are denned in Ref. |5T| . 

We have computed the functions G c (z) and S E (z) using the numerical technique pre- 
sented in Ref. ||38|| . In the resummation we have used the seven-loop results of Ref. |57] that 
allow the extension of the series expansions of j(g) and k(g) by one order. If 

3(9) = J2^ n i = E ^ ( A102 ) 

n n 



we derive from Ref. 57 



j7 = 0.0996888, k 7 = -0.00190671. (A103) 

The results for G c (z) are well fitted by 

G c (z) = (1 + 50.79365z + 508.54282 2 + 5929.475z 3 + 10937.032 4 ) - 07875 . (A104) 

For z — > this expression gives G c (z) ~ l+4z, in agreement with the perturbative expansion 
(|S8|), while for z — > 00 we have 

G c (z) = 2.08 z 2 ^ 2 (1 + 0.089* -1 + 0(^ 2 )). (A105) 

This expansion is in reasonable agreement with Eq. ( |A93j ), keeping into account that A w 
1/2. However, while the leading term is close to the estimate ( |A95| ), the correction differs 
significantly from Eq. ( A97 ). A simpler expression, that is however more accurate in the 
Wilson-Fisher region, is 

G c (z) = (1 + 38.0952z + 276.8442 2 + 1073.172 3 ) - 105 . (A106) 

For z — > 00 it behaves as 

G c (z) w 2.08 z 2 ^ 2 (1 + 0.027Z" 1 ), (A107) 
in agreement with the asymptotic expansion (|A93j ) and the numerical values ([A95Q , ( |A97|) . 



For Se(z) we find that the expression reported in Ref. |pl 



S E (z) = (1 + 7.61182 + 12.05135^) ai/5it>b , (A108) 

fits the data extremely well and correctly reproduces the asymptotic behavior for z — > and 

z — >• 00. 

We mention that a different interpolation appears in Ref. [ |55[] , based on the five-loop 
computations of Schloms and Dohm p9| . Setting z = 6.95? in the formulae of Ref. |89[ 
order to reproduce the correct behavior for z — > 0, we obtain 
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S E ( Z ) = (1 + 7.4074^ + 10.913^) aiS , (A109) 

in good agreement with Eq. ( |A108j ). Finally, from the results reported in Ref. we can 
estimate the ratio g c \j s E \- We find g c \j s E \ ~ 0.26, in reasonable agreement with Eq. (|A99|) . 
Let us finally notice that if we use the estimate ( |A29[ ) we would obtain g c \j s E \ ~ 0.23. 
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TABLE II. Monte Carlo results. Here /3 mf = (V R - l) -1 . 



p 


n 


log 


log (Cn/3 ,") 
' mT ' 


2 


20 


4.094106(9) 


-1.570800(12) 




30 


4.556891(21) 


-2.52151(4) 




40 


4.886886(17) 


-3.48884(3) 




80 


5.68682(3) 


-7.42459(6) 




160 


6.49153(6) 


-15.39944(14) 




320 


7.29966(12) 


-31.45439(30) 




640 


8.10973(24) 


-63.67150(64) 




1280 


8.92194(54) 


-128.2124(16) 


3 


120 


6.50585(3) 


-5.20178(6) 




160 


6.82940(15) 


-7.06537(30) 




240 


7.28750(5) 


-10.81482(13) 




320 


7.61466(26) 


-14.58095(62) 




480 


8.07689(10) 


-22.13565(27) 




960 


8.87276(19) 


-44.87624(56) 




1920 


9.67317(37) 


-90.4588(12) 




3840 


10.47830(81) 


-181.7297(31) 


4 


1200 


9.39403(32) 


-30.4356(12) 




1280 


9.46739(15) 


-32.4988(5) 




1920 


9.92722(12) 


-49.0221(4) 




2560 


10.25489(28) 


-65.5637(10) 




3840 


10.71901(22) 


-98.6670(8) 




5120 


11.04855(54) 


-131.7917(21) 




7680 


11.51585(45) 


-198.0564(17) 




10240 


11.8470(11) 


-264.3489(44) 




15360 


12.31571(93) 


-396.9388(36) 




30720 


13.1234(22) 


-794.8075(75) 


5 


600 


8.88068(10) 


-8.81049(25) 




1040 


9.484413(43) 


-15.51334(13) 




1200 


9.64205(17) 


-17.95699(50) 




2080 


10.25375(8) 


-31.42686(26) 




2400 


10.41268(28) 


-36.3313(10) 




4160 


11.03218(14) 


-63.3411(6) 




8320 


11.81856(25) 


-127.2624(11) 




16640 


12.61152(54) 


-255.2014(24) 




33280 


13.4103(11) 


-511.1808(48) 




66560 


14.2166(31) 


-1023.251(12) 


6 


500 


8.93513(4) 


-4.61781(8) 




800 


9.43841(8) 


-7.52107(20) 




1000 


9.67871(6) 


-9.46419(16) 




1600 


10.18886(14) 


-15.31451(40) 




2000 


10.43242(10) 


-19.22308(33) 




3200 


10.94927(23) 


-30.97376(82) 




4000 


11.19612(18) 


-38.81610(66) 




8000 


11.96920(31) 


-78.0846(13) 




16000 


12.75136(54) 


-156.7110(27) 




32000 


13.5400(11) 


-314.0585(56) 
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TABLE III. Monte Carlo results. Here /3 mf = (V R - l)" 1 . 



p 


n 


log #2 


log(c n /3-™) 


7 


100 


7.500434(17) 


-0.543304(13) 




200 


8.214646(25) 


-1.164668(26) 




400 


8.93561(4) 


-2.44010(5) 




800 


9.66443(6) 


-5.03235(11) 




1600 


10.40211(9) 


-10.26761(23) 




3200 


11.14984(16) 


-20.79866(51) 




6400 


11.90757(25) 


-41.9304(10) 




12800 


12.67547(43) 


-84.2725(21) 




25600 


13.45153(68) 


-169.0420(42) 




51200 


14.2366(14) 


-338.6733(87) 


8 


750 


9.80702(3) 


-3.29563(5) 




1500 


10.53428(5) 


-6.74019(10) 




3000 


11.27031(7) 


-13.67826(20) 




6000 


12.01591(12) 


-27.61282(42) 




12000 


12.77156(19) 


-55.55008(85) 




24000 


13.53685(34) 


-111.5009(17) 




48000 


14.31143(59) 


-223.4866(35) 


9 


1040 


10.34095(12) 


-3.35462(21) 




2080 


11.06434(17) 


-6.8409(4) 




4160 


11.79586(27) 


-13.8574(9) 




8320 


12.5369(4) 


-27.9440(18) 




16640 


13.2868(6) 


-56.1807(36) 




33280 


14.0491(10) 


-112.7270(71) 




66560 


14.8188(20) 


-225.901(14) 


10 


1500 


10.89765(8) 


-3.66081(13) 




2000 


11.19598(5) 


-4.91804(11) 




3000 


11.61839(11) 


-7.44241(26) 




4000 


11.91986(8) 


-9.97288(22) 




6000 


12.34684(17) 


-15.04661(52) 




8000 


12.65209(13) 


-20.12732(44) 




16000 


13.39353(20) 


-40.49032(91) 




32000 


14.14511(32) 


-81.2802(18) 




64000 


14.90648(72) 


-162.9324(51) 


12 


1500 


11.21510(8) 


-2.19646(10) 




2000 


11.50974(10) 


-2.95223(16) 




3000 


11.92649(14) 


-4.47024(21) 




4000 


12.22324(16) 


-5.9926(4) 




6000 


12.64352(20) 


-9.0463(4) 




8000 


12.94315(25) 


-12.1053(8) 




12000 


13.36771(30) 


-18.2350(9) 




16000 


13.67113(37) 


-24.3704(16) 




32000 


14.40765(57) 


-48.9504(31) 
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TABLE IV. Determination of (3 Cy R for several values of p. The reported results are obtained 
by fitting the numerical data with n > n m \ n . Fit (a): VR/3 e fT,il( n ) = Vr(3 c ,r + a/n + W? 3 /n L5 . Fit 
(b): V R (3 cS:R (n) = V R (3 C , R + a/n + bR 3 /n 15 + cR 6 /n 2 . 



p 


^min 


VrPcR 


a 


b 


c 


Fit (a) 


2 


40 


1.152388(1) 


0.0725(4) 


0.1044(16) 




3 


80 


1.0656569(3) 


0.0142(3) 


0.0360(5) 




4 


150 


1.0342660(1) 


-0.0040(3) 


0.0174(3) 




5 


260 


1.0199238(1) 


-0.0111(1) 


0.0100(1) 




6 


400 


1.0125749(1) 


-0.0138(4) 


0.0066(2) 




7 


800 


1.00840580(4) 


-0.0164(3) 


0.0053(1) 




8 


1500 


1.00588845(3) 


-0.0174(4) 


0.0044(1) 




9 


1040 


1.00427552(9) 


-0.0140(10) 


0.0026(2) 




10 


3000 


1.00320022(5) 


-0.0167(13) 


0.0031(3) 




12 


2000 


1.00192307(7) 


-0.0080(11) 


0.0010(1) 




Fit (b) 


2 


10 


1.152388(1) 


0.0711(1) 


0.1211(3) 


-0.01139(6) 


3 


15 


1.0656574(3) 


0.0121(1) 


0.0442(1) 


-0.00182(1) 


4 


40 


1.0342661(1) 


-0.0057(1) 


0.0217(1) 


-0.000589(6) 


5 


70 


1.0199239(1) 


-0.0130(1) 


0.0132(1) 


-0.000274(3) 


6 


100 


1.0125750(1) 


-0.0158(4) 


0.0090(2) 


-0.000152(6) 


7 


140 


1.00840566(10) 


-0.0159(4) 


0.0062(1) 


-0.000081(3) 


8 


300 


1.00588839(5) 


-0.0176(5) 


0.0055(2) 


-0.000071(3) 


9 


250 


1.00427551(10) 


-0.0153(10) 


0.0037(2) 


-0.000037(3) 


10 


450 


1.00320011(4) 


-0.0143(5) 


0.0030(1) 


-0.000028(1) 


12 


400 


1.00192275(6) 


-0.0028(5) 


0.0006(1) 


-0.000003(1) 
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FIG. 1. Estimates of T2, e fr versus 1/R. The reported points correspond to 
p = 3, 4, 5, 6, 7, 8, 10, 12. The line is the best fit: T 2 , e ff = -0.00814 + 0.00277/i?. The errors 
on the data take into account the error on f3 Cj R and on the constant t±. 
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FIG. 2. Estimates of log c rai j? versus n = nR~ 6 . The reported points correspond to 
p = 3,4, 5, 6, 7, 8, 10, 12. The line is the field-theoretical prediction. In the upper graph we re- 
port all points, in the lower one only those satisfying n > Vr/2. 
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FIG. 3. Estimates of c U) r/ g c ,th{n) versus n = nR~ 6 . The reported points correspond to 
- 3, 4, 5, 6, 7, 8, 10, 12. The errors only take into account the error on c nt R. 
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FIG. 4. Estimates of A c;n ^ versus n = nR~ 6 . The reported points correspond to 
: 3,4,5,6,7,8,10,12. 
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OR" 6 

FIG. 5. Estimates of Q n ,R (upper graph) and of Aq- h ^ (lower graph) versus n = nR~ 6 . The 
reported points correspond to p = 3, 4, 5, 6, 7, 8, 10, 12. The continuous line in the upper graph is 
the theoretical prediction. 
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FIG. 6. Estimates of 7 e ff(n, R) versus n = ni?~ 6 . The reported points correspond to 
p = 2, 3, 4, 5, 6, 7, 8, 10, 12. The continuous line is the theoretical prediction. 
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FIG. 7. Estimates of E^ R /(6n) (upper graph) and of /S.E- n ,R (lower graph) versus n = nR~ e . 
The reported points correspond to p = 3, 4, 5, 6, 7, 8, 10, 12. The continuous line in the upper graph 
is the theoretical prediction. 
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FIG. 8. Estimates of u e g(n,R) versus n = nR~ 6 . The reported points correspond to 
p = 2, 3, 4, 5, 6, 7, 8, 10, 12. The continuous line is the theoretical prediction. 
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FIG. 9. The effective exponents 7 e fj (upper curve) and u e g (lower curve) n = nR 6 , obtained 
from the phenomenological expressions (|102[), (|103|), for several values of p. 
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